Planar defects and the fate of the Bragg glass phase of type-II superconductors 
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It is shown that the Bragg glass phase can become unstable with respect to planar defects. A 
single defect plane that is oriented parallel to the magnetic field as well as to one of the main axis 
of the Abrikosov flux line lattice is always relevant, whereas we argue that a plane with higher 
Miller index is irrelevant, even at large defect potentials. A finite density of parallel defects with 
random separations can be relevant even for larger Miller indices. Defects that are aligned with 
the applied field restore locally the flux density oscillations which decay algebraically with distance 
from the defect. The current voltage relation is changed to lnV(J) ~ —J -1 . The theory exhibits 
some similarities to the physics of Luttinger liquids with impurities. 

PACS numbers: 74.25.Qt, 61.72.Mm, 72.15.Rn 



Typc-II superconductors have to contain a certain 
amount of disorder to sustain superconductivity: the dis- 
order pins magnetic flux lines, hence preventing dissipa- 
tion due to their motion |lj, |2|. For some time it was 
believed that disorder due to randomly distributed impu- 
rities destroys the long range translational order (LRO) 
of the Abrikosov flux line lattice |3( ■ More recently it was 
shown that the effect of impurities is much weaker result- 
ingin a phase with quasi-LRO, the so-called Bragg glass 
Lit Ll lil til I n tn i s phase the averaged flux line density is 
constant but a remnant of its periodic order is seen in the 
correlations of the oscillating part of the density which 
decay as a power law. Experimental signatures of this 
phase have been observed [8j. An important feature of 
the Bragg glass is the highly non-linear current-voltage 
relation related to the flux creep which is of the form 
InV(J) ~ — J -1 / 2 so that the linear resistance vanishes. 

Although much of the original transport data on flux 
creep in high-T c superconductors was discussed in terms 
of point disorder (see e.g. |9() it was realized later that 
many samples included planar defects like twin planes or 
gra in boundaries which masked the Bragg glass behavior 
jinl | . Indeed, in clean samples planar defects lead to much 
more pronounced pinning phenomena than point disorder 
because of stronger spatial correlations |2,lH|- However, 
the generic experimental situation is a mixture of point 
disorder and planar defects, a case which has not been 
studied theoretic ally in the context of the Bragg glass 
yet; see however [l2lll3j . 

It is the aim of the present paper to consider exactly 
this case, i.e. the question of the influence of planar de- 
fects in the Bragg glass phase. Our key results are as fol- 
lows: a necessary condition for a planar defect to become 
a relevant perturbation is that it is oriented parallel to 
the magnetic flux. In this case its influence on the Bragg 
glass phase can be characterized by the value of a single 
parameter g = §77 (a / '£) 2 which depends both on the ex- 
ponent 77 describing the decay of the density correlations 



in the Bragg glass phase and on the orientation of the 
defect with respect to the flux line lattice, a and I are 
the mean spacing of the flux lines in the bulk and the 
distance between lattice planes of the Abrikosov lattice 
parallel to the defect, respectively. The defect is relevant 
for an infinitesimal weak defect potential if g < 1 which 
is realized if and only if the defect plane is parallel to 
one of the main crystallographic planes of the flux line 
lattice. In the vicinity of the (relevant) defect the density 
profile shows periodic order with an amplitude decaying 
algebraically (with exponent g) with the normal distance 
from the defect. The current voltage relation for volt- 
age drops perpendicular to a defect plane is of the form 
In Vd(J) ~ — J -1 . Correlations are destroyed across (rel- 
evant) defects. For g > 1 a weak (and presumably even 
a strong) defect is irrelevant and the density profile de- 
cays with a larger exponent 2g — 1 > 1 . All defects which 
are tilted against the magnetic flux are irrelevant and the 
flux density oscillations decay exponentially. For a finite 
density of parallel defect planes (with random distances) 
the Bragg glass is destroyed for g < g c with g c > 3/2. 

It is worthwhile to mention that some aspects of 
our results are similar to other theories at their crit- 
ical dimension like 2-dimensional classical or (1+1)- 
dimensional quantum models. Adding a planar defect 
in the Bragg glass resembles the presence of a line defect 
in the classical or a frozen impurity in the quantum case 
|l4J lla . llfjl |l?J when the coupling constant g is identi- 
fied with temperature or the Luttinger liquid parameter, 
respectively. The periodic order seen in the vicinity of 
the defect plane has its counterpart in the Friedel os- 
cillations in Luttinger liquids close to an isolated impu- 
rity. Whereas in those cases the relevance of a defect 
(i.e. an impurity) can be changed by tuning the inter- 
action strength or the temperature, respectively, in the 
present case a change of g can be accomplished by chang- 
ing the orientation of the defect with respect to the flux 
line lattice. The current voltage relations are however 



different from the 2-dimensional cases. Experimentally 
micro-twinned crystals with one direction of twin planes 
have been studied, e.g., in Ref. [lSj . 

First, we consider the effect of a single planar defect 
in a system of interacting flux lines which are pinned by 
randomly distributed impurities. Since we are interested 
in the behavior on large length scales, it is appropriate 
to describe the interacting flux lines in terms of a con- 
tinuum elastic approximation with a displacement field 
u(r). The Hamiltonian H = Hq + Hd includes then the 
clastic energy and the impurity interaction, 

H a = - /"dzd 2 x{ Cll (V x -u) 2 +c 66 (V x xu) 2 + 

+ c 44 (d z uf + £ , p in (u,r)| , (1) 

and the coupling energy Hd to the defect, see be- 
low. The pinning energy can be written as E pin — 
U(r)p(r,u), where U(r) is a potential arising from ran- 
domly distributed impurities and r = (x, z). The lo- 
cal flux line density can be expressed as p(r, u) = 
p {-V x u(r)+£ G elG[x ~ u(r)1 } 0. Here Po = 2/V3a 2 , 
and G = G TOn = mbi + nbj is a vector of the reciprocal 
lattice with integer m, n 23]. The set of the six smallest 
reciprocal lattice vectors will be denoted by {Go}. 

The energy of a planar defect has the same form as 
-Epin but with U(r) replaced by — vS(r ■ ud — 5) where v 
is the strength of the defect potential, ri£> and S denote 
the normal vector of the defect plane and its distance 
(along no) from the origin of the coordinate system, re- 
spectively. This gives 



H D = vpo fd 3 r <5(r • n D - S) { V x u(r) - ^ < 

"" G#0 



iG[x-u(r)]l 

(2) 

Without Hd this model has been studied in detail 
using different approaches [J, la, |g, L|- It was shown 
that thermal fluctuations are irrelevant (zero tempera- 
ture fixed point) and that the pinned flux line lattice 
exhibits a power law decay of the translational order pa- 
rameter V[ / g(x, z) = e " lGu ( x ' z ') t similar to pure 2D crys- 
tals at finite temperatures. In particular, ("Jg) = but 
(* G (x,0)*_ G (O,0)) - Ixl"^, where (...) denotes both 
the thermal and disorder average and tjg = r](G/Go) 2 - 
From a Gaussian variational treatment in d — 3 dimen- 
sions follows t] = 1 |6j whereas a functional renormal- 
ization group analysis in d = 4 — e dimensions yields 
a non-universal exponent r\ that varies with the elas- 
tic constants of the vortex lattice pj. Extrapolating to 
d = 3 [24|, one finds only a very weak variation with 
1.143 < ?/ < 1.159 0. Since (despite of the glassy nature 
of the phase) the structure factor shows still Bragg peaks 
the notation Bragg-glass was coined |o] ■ Next we discuss 
the influence of Hd ■ In order to integrate over the delta 
function in Eq. Pp. it is convenient to introduce an ex- 
plicit parameterization for the position vector vd of the 




FIG. 1: Triangular flux line lattice with vectors of the direct 
(ai and &2) and the reciprocal lattice (bi and b2) and the 
orientation of the defect plane. 

defect plane which obeys r^i • n^ = 5. With 

r D = (x D ,z D ) + 5n D , z D = tcos/3 (3) 

X£> = (ssina — t cos a sin /3, scosa + t sinasin/3) 

we introduce in-plane coordinates s, t, and two angles 
a and (3 which determine the rotation of the plane with 
respect to the y- and z-axis, respectively (see Figure HJ. 
The defect energy reads then 



H D = vp / dtds {V x u(r D ) - 



G^O 



iG[i5n D +x D — u(r_o)] "I 

(4) 

Since the displacement field u(ro) varies slowly on the 
scale of the flux line lattice constant a, the integrals 
over s and t vanish for all G except those for which 
the oscillatory factor e lGxD becomes one (for all s, t). 
This condition can be satisfied only if sin/3 — 0, i.e., 
if the defect plane is parallel to the applied magnetic 
field. There remains a second condition for the angle 
a which results from the constraint that G has to be 
perpendicular to xd- Expressing the defect plane (for 
sin/3 = 0) as X£> = (ciai — C2Si2)s/a where a.;bj = 27r<5.y, 
this results in the condition m/n = C2Jc\. Hence if 
C1/C2 is irrational the effect of the defect plane is al- 
ways averaged to zero. On the other hand, for ratio- 
nal C2/C1 we may choose rriD,nD to be the smallest co- 
prime pair with C2/C1 = m.DJnD- Then mo, riD are 
the Miller indices of the defect plane and only those G 
which are integer multiples of Gd = G„ lE ,„ D contribute 
in Eq. Q. In the following, we will concentrate on the 
contribution from these G-vectors only. The flux line 
lattice planes (of the ideal lattice) parallel to a defect 
plane with Miller indices to.d, ud have a separation of 
£ = ^a/ \/m D + uidud + n D and hence Gd = 2tt/£- 

Using the result for the average of ^G(r) in the Bragg 
glass phase [?| , one finds for the disorder averaged defect 



energy Ho on scale L 

(Wd)o 



^2v k cos(kSG D ) ( — 
fc=i ^ a 



2-k 2 g 






where (. . .)o denotes the average with Ho- The linear 
gradient term in J3J vanishes after averaging [25j. The 
coefficients cos(kSGD) reflect the periodicity of the de- 
fect energy under translations by £ normal to the plane. 
The Gaussian approximation used here is believed to be 
correct to order e |l9| . It is important to remark that the 
i-dcpcndence of Eq. (JSJ) holds only on length scales larger 
than the positional correlation length L a w Lg(a/£) 1/ ^ rm 
where £ is the maximum of the coherence and the disor- 
der correlation length and L^ denotes the Larkin length 
on which the typical flux line displacement is of the or- 
der £. £ rTO w 0.175 is the roughness exponent of the 
elastic distortions on scales smaller than L a (in the so- 
called random manifold regime). There the correlations 
of V^gOO decay as a stretched exponential and the ef- 
fect of the defect plane is reduced by disorder fluctua- 
tions. Due to these fluctuations on intermediate length 
scales the initial value of the defect strength is reduced 
to Vk ~ v{L a /a) 2 e~ c ^ GDka ^ where c is a constant. 

To linear order in Vk, the RG flow equation of the Vk 
is obtained by comparison of the defect energy scaling in 
in Eq. (JjjJ with the scaling of Hq at the Bragg glass fixed 
point, yielding 



diik/dlnL = (1 — k g) v k 
Hence v\ is a relevant perturbation for g 



(6) 



rj{m 2 D +m D n D + n 2 D ) < 2 or 



g < 1, i.e., if 



0.66 a, (7) 



which is compatible only with £ — V3a/2 « 0.87a. Hence 
the defect plane must be oriented parallel to one of the 
three main crystallographic planes of the flux line lattice 
(i.e. cos 2/3 = cos 6a = 1). 

The transition described by Eq. ijfjj occurs not in the 
bulk but on the defect plane. Hence one can develop 
an effective theory on the defect which could be used to 
describe also stronger defect potentials. Since the de- 
fect couples only to the normal displacement u±(td) = 
ri£>u(r£>) ° n the defect plane, we would like to integrate 
out uj_ outside the defect and u — u±no across the entire 
sample. This integration is facilitated by the reasonable 
assumption of an effective Gaussian theory for the defect- 
free system at the Bragg glass fixed point 0. We find 
that u± on the defect has long-ranged elasticity and is 
described by the effective Hamiltonian (compare |16l [l7j 
for a corresponding procedure in the clean case) 



n 



2D 



f /^qlqll^l 2 



(8) 



U 2 r D 



— cos(^ - a) + — cos(<p) >■ 



where ip(ro) = 2iru±(ro)/£ and q is the in-plane mo- 
mentum, a is a random phase which is uncorrelated and 
uniformly distributed. The amplitude of the random po- 
tential has been chosen here as to reproduce the proper 
Gaussian replica theory of the Bragg glass for vi = 0. 
The elastic constant K depends on the bulk elastic mod- 
uli, the angle a and the direction of q. The model of 
Eq. (JHJ shows a transition at g = g c (vi) with g c (0) = 1 
in agreement with our previous considerations. In the 
present case g = r)(m D + mono + n D )/2 can only be 
changed in finite steps by changing the orientation of the 
defect plane. Thus at small v\ only the defects paral- 
lel to the three main crystallographic planes are relevant 
(with m D + mono + n D = 1 ). Due to the long-ranged 
elasticity we expect from our analysis of a similar model 
|20| that even at large V\ one has g c {vx) = 1. 

Next we study the order of the flux lines in the vicinity 
of the defect. Since the superconducting order is reduced 
in the defect plane it is plausible to assume that v > 0. 
The potential of a relevant defect growth under renormal- 
ization and effectively decouples the two half-spaces. On 
large scales, we can then impose Dirichlct boundary con- 
ditions for the normal flux line displacement u± on the 
defect plane. This allows us to determine the boundary 
induced modifications of the flux density via the method 
of images. We find that (uj_(r)) = i([u ± (r)-u ± (r m )] 2 ) , 
where r m is the mirror image of r with respect to the de- 
fect plane. From this we obtain immediately that with 
increasing distance L± = |noX — S\ from the defect, the 
slowest oscillations of the flux line density decay for a 
relevant defect plane (g < 1) as 



(p(x,z,u)) - ( j^- 



cos{G (L ± ±S)). 



(9) 



Hence, a single relevant defect plane yields a long-ranged 
restoration of the order parameter. The oscillations of 
the density resemble Friedel oscillations observed in Lut- 
tinger liquids close to an isolated impurity |17|. This 
similarity is substantiated by considering the decay of 
the density oscillations if the defect plane is irrelevant, 
i.e., for p > 1. Then the defect potential decreases un- 
der renormalization and lowest order perturbation the- 
ory can be applied. Such an approach takes into account 
that the defect strength decreases as v(L) ~ L 1 ^ 9 and 
hence we obtain Eq. © with g replaced by 2g — 1 > 1 , 
in close analogy to the (1 + l)-dimcnsional counterpart 

0EG1. 

If the defect plane is not parallel to the applied mag- 
netic field, i.e. sin/3 ^ 0, then the irrelevance of the 
defect should be also observable in the Friedel oscilla- 
tions of the flux line density which should become then 
short-ranged. If g > 1 we can use perturbation theory to 
obtain the density oscillations since the defect strength 
deceases for large length scales, i.e., for sufficiently small 
sin/3. From this analysis we find for the slowest oscilla- 



tions normal to the defect (with 5 = 0) 



(p(x,z,u)) 



r,2 



-L±/i 



cos (GdL± cos/3) 



(10) 

for large separations L± ^> £ with a characteristic length 
scale £ = 1/(G.d| sin/3|). Hence, if the tilt angle (3 ap- 
proaches zero, £ diverges, and the Friedcl oscillations for 
L± <C £ are described by Eq. (|10|) with £ replaced by 
Lj_, hence in agreement with our findings for (3 = 0. Al- 
though our perturbative analysis is restricted to g > 1 we 
expect that for any value of g the Friedel oscillations de- 
cay exponentially beyond distances of order £. Since the 
renormalization of the defect potential is cutoff on the 
scale £, even for a relevant defect perturbation theory is 
expected to be justified if both the bare defect potential 
is sufficiently weak and sin/3 is not too small. 

To find the current voltage relation for a relevant de- 
fect we consider the most interesting case where the cur- 
rent J is parallel to the defect and normal to the mag- 
netic field B. Flux creep in the presence of pinning forces 
arises via formation of critical droplets [2| ■ In the present 
case of a defect plane the droplet is characterized by 
Gdu^{td) = and 2-k outside and inside of the droplet, 
respectively. Since on scale L a typical distortions in the 
Bragg glass phase are of the order a, the droplet volume 
is ~ L 2 L a and the volume energy gain in the droplet is 
of the order ~ —JBL 2 L a which has to be compared with 
the the surface energy loss of the order ~ v x ' ~ 9 ' LL a 
where we have included the renormalization of v\ . From 
the balance of the two terms follows a critical droplet 



size L„ 



,i/2(i-s) 



(L a JB) 1 . Droplets of this size have 
a free energy ~ v x (JB)^ 1 which determines via the 

Arrhenius law the creep velocity due to thermal activa- 
tion, yielding the voltage drop normal to the defect plane 



F~exp[-(C /(TJ))], Co 



,1/(1-2) 



zr 



(ii) 



Thus the flux creep across the defect is much slower than 
in the Bragg glass phase. If for weak defects g — > 1 — , Co 
becomes small, hence reducing the applicability of this 
formula to extremely small currents. 

Crystals usually contain either two orthogonal families 
( "colonies" ) of parallel twin boundaries or a single family 
of parallel twin planes. We focus on the latter case, which 
for impure s amp les has been studied so far only in (1+1) 
dimensions |21|. We consider the typical situation of 
planes with random distances Sj that are aligned with the 
applied field. Each plane contributes Hd(Sj) of Eq. (@J 
with Sj replacing S. For the averages over the Sj one finds 

{LjHDiSj)) = but <E,W?>(^)) 1/2 ~ L b/2 ~ 9 from the 
term with k = 1. Since TLq ~ L, a family of weak de- 
fect planes is relevant and destroys the Bragg glass for 
g < g c = 3/2, i.e., the planes must be oriented with a 



main plane of the flux lattice. Contrary to a single de- 
fect plane, we expect that a larger defect strength yields 
an increased g c > 3/2 rendering additional defect orien- 
tations relevant. However, a description of the transition 
at strong coupling and the localized fixed point describ- 
ing relevant defect planes is not available at present [22| ■ 
Only for defect planes of equal distance we expect long 
range order in the direction perpendicular to the planes. 
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A direct computation for d = 3 yields a value for g that 
depends not only on cn/cee but on the two ratios ch/cm 
and C66/C44. However, the range of variation of g remains 
unchanged to order e. 

The coupling of V x u(r) to the defect potential for a de- 
fect parallel to the z-axis can be eliminated by the trans- 
formation u — > u + V x </>(x) with </>(x) a scalar field which 
obeys V^ = ^°-S(x-iid — 5). This transformation does 
not change the terms ~ egg, C44 in Eq. Q but shifts the 
flux line density on the defect. 



